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How robust are predictions of galaxy clustering? 
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1 INTRODUCTION 



Galaxy formation is an inefficient process, with only a few per- 
cent of the available baryons in t he Universe found in a "con- 
densed" state as stars or cold gas dBalogh et all 1200 it ICole et alj 
l200ll ; iMcCarthv et al.ll2007l) . The fraction of condensed baryons 
also varies strongly with halo mass, as a result of the interplay be- 
tween the processes which participate in galaxy formation, such as 
gas cooling, star formation, heating of the interstellar medium by 
super novae and the impact on the host galaxy of energy released by 
AGN teaughll200dlBensonll2oToh . Physical calculations of galaxy 
formation attempt to model all of these processes to predict the 
number of galaxies hosted by a dark matter halo and their proper- 
ties. The aim of this paper is to assess how robustly current models 
can predict the number of galaxies which are hosted by dark matter 
haloes of different mass. By focusing on a subset of the predictions 
possible with current models and through selecting galaxies based 
on their intrinsic properties rather using more direct observational 



ABSTRACT 

We use the Millennium Simulation database to compare how different versions of the Durham 
and Munich semi-analytical galaxy formation models populate dark matter haloes with galax- 
ies. The models follow the same physical processes but differ in how these are implemented. 
All of the models we consider use the Millennium N-body Simulation; however, the Durham 
and Munich groups use independent algorithms to construct halo merger histories from the 
simulation output. We compare the predicted halo occupation distributions (HODs) and corre- 
lation functions for galaxy samples defined by stellar mass, cold gas mass and star formation 
rate. The model predictions for the HOD are remarkably similar for samples ranked by stellar 
mass. The predicted bias averaged over pair separations in the range 5 - 25/?~'Mpc is con- 
sistent between models to within 10%. At small pair separations there is a clear difference in 
the predicted clustering. This arises because the Durham models allow some satellite galaxies 
to merge with the central galaxy in a halo when they are still associated with resolved dark 
matter subhaloes. The agreement between the models is less good for samples defined by cold 
gas mass or star formation rate, with the spread in predicted galaxy bias reaching 20% and 
the small scale clustering differing by an order of magnitude, reflecting the uncertainty in the 
modelling of star formation. The model predictions in these cases are nevertheless qualita- 
tively similar, with a markedly shallower slope for the correlation function than is found for 
stellar mass selected samples and with the HOD displaying an asymmetric peak for central 
galaxies. We provide illustrative parametric fits to the HODs predicted by the models. Our 
results reveal the current limitations on how well we can predict galaxy bias in a fixed cos- 
mology, which has implications for the interpretation of constraints on the physics of galaxy 
formation from galaxy clustering measurements and the ability of future galaxy surveys to 
measure dark energy. 
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criteria, we will be able to make a cleaner comparison between the 
models. 

For this comparison we use publicly available galaxy cata- 
logues produced by two groups who have independently developed 
semi-analytical models of galaxy formation. Semi-analytical mod- 
els attempt to calculate the fate of the baryonic component of the 
universe, in the context of the hierarchical growth of structure in the 
dark matter. These models use differential equations to describe the 
processes listed in the opening paragraph. Often, these processes 
are poorly understood and nonlinear, so the equations contain pa- 
rameters. The values of the parameters are set by comparing the 
model predictions to a selection of observational data, and adjust- 
ing the parameter values until an acceptable match is obtained. Cur- 
rently, and for the foreseeable future, semi-analytical modelling is 
the only technique which can feasibly be used to populate large 
cosmological volumes of dark matter haloes with galaxies to ob- 
tain predictions for galaxy clustering out to scales of several tens of 
megaparsecs. 

In this comparison we use galaxy formation models 
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which have been ru n using the Millennium N-body simulation 
dSpringel et al]|2005t) . We consider a range of models run by the 
"Durham" and "Munich" groups (listed in the next section) using 
the outputs from the dark matter only Millennium Simulation. The 
two groups have their own algorithms for constructing merger his- 
tories for dark matter haloes and different implementations of the 
physics of galaxy formation. Since the models populate the same 
dark matter distribution with galaxies this provides an opportunity 
to look for any systematic differences in the predictions for the 
galaxy content of dark matter haloes. 

The conclusions of this comparison will tell us how robust the 
predictions of current models are, given the uncertainties in the un- 
derlying physics. This is important for assessing how useful mea- 
surements of galaxy clustering are for constraining the physics of 
galaxy formation. If the models purport to follow the same pro- 
cesses, yet predict different numbers of galaxies per halo, then this 
limits what we can leam from clustering at present. As well as im- 
proving our understanding of physics, modelling galaxy clustering 
and how it relates to the clustering of the underlying dark matter, 
called galaxy bias, is also impo rtant for future galaxy surveys which 
aim t o measure dark energy dLaureiis et all boilt [Schlegel et al] 



Galaxy bias is a systematic which limits the performance 
of large-scale structure probes of dark energy. If we can model bias 
accurately, then this systematic can be marginalized over. 

This paper is organized as follows. In Section 2 we give a brief 
overview of semi-analytical modelling and state which models we 
are comparing. Section 3 describes some preparatory work for the 
comparison, which involves post-processing of the catalogues and 
describes how we construct our samples. The main results of the 
paper are in Section 4 and our conclusions are presented in Section 
5. The Appendix describes parametric fits to the halo occupation 
distribution predicted by the models. 



2 THE GALAXY FORMATION MODELS 

We compare the predictions of five different semi-analytical mod- 
els of galaxy formation which are publicly available from the Mil- 
lennium Archive in GarchinjQ and its mirror in Durham^]. The 
models are all set in the cosmological context of the Millennium 
N-body simulation o f the hierarchical cl ustering of matter in a 
A-CDM cosmology dSpringel et alj[2005h . The models were pro- 
duced by two independent groups of researchers, and correspond 
to "best bet" models released to the community since 2006. Two 
of the mo dels are generally r eferred to under the label of "Durham 
models" dBower et al]|2006r . iFont et al.ll2008h and will be referred 
to in plots, respectively, as "Bower 20 06" and "Font 2008", and 
the other three are "Munich" model s dDe Lucia & Blaizotj |2007l ; 
iBertone et afll2007l; iGuo et alj [2oTlT) . which will be labelled as 
"DeLucia2007", "Bertone2007" and "Guo2011", respectively. 

The models all aim to follow the main physical processes 
which are believed to be responsible for shaping the formation and 
evolution of galaxies. These are: (i) the collapse and merging of 
dark matter haloes; (ii) the shock heating and radiative cooling of 
gas inside dark matter haloes, leading to the formation of galaxy 
discs; (iii) quiescent star formation in galaxy discs; (iv) feedback 
from supernovae (SNe), accretion of mass onto supermassive black 
holes and from photoionization heating of the intergalactic medium 



(IGM); (v) chemical enrichment of the stars and gas; (vi) galaxy 
mergers driven by dynamical friction within common dark matter 
haloes, leading to the formation of stellar spheroids, which may 
also trigger bursts of star formation. However, the implementations 
of all of these processes differ between the models. These differ- 
ences even include the first step listed above of the construction of 
halo merger trees from the N-body simulation. The modelling of 
the above processes is uncertain and the resulting equations often 
require parameter values to be specified. The models differ in how 
these parameters are set, as the different groups assign different im- 
portance to reproducing various observational datasets. 

It is not our intention to present a comprehensive descrip- 
tion of the models and their differences. Full details of the mod- 
els can be found in the references given above and in earlier pa- 
pe rs by each group. T he Durham model, GALFORM, was introduced 
by ICole et al] d2000h a nd extended, for the purposes of the mod- 
els considered here by iBenson et alj d2003l). The Mu nich model, 
LGALAXIES, was introdu ced by ISpringel et al l J200lh and devel- 
oped in a series of pa pers dDe Lucia et al .l2004l;ICroton et al.l2006l ; 
|Pe Lucia et al] 20061) . 

Instead, as we present our results and try to interpret the level 
of agreement between the model predictions, we will discuss var- 
ious components of the models which we believe are responsible 
for any differences. 



3 PRELIMINARIES : PREPARATION FOR A 
COMPARISON BETWEEN MODELS 

Having downloaded the galaxy catalogues from the Millennium 
Archival, in order to carry out a robust and meaningful comparison 
between the model predictions, it is important to account for any 
differences in definitions of properties and to set up well defined 
samples. 

There are two aspects we need to homogenize for our compar- 
ison: the definition of mass used to label dark matter haloes and the 
values of galaxy properties used to define samples. We deal with 
each of these in turn below. We close this section by discussing a 
relabelling of some of the halo masses that we found necessary in 
the Munich models, due to differences in the algorithms used to 
construct the halo merger histories. 



3.1 Definition of halo mass 

The Durham models list halo masses (archive property name 
"mhalo") derived fro m the "DHalos" me rger tree construction, 
which is described in Merson et alj d2012l) (see also Jiang et al. 
in preparation). The algorithm is designed to ensure that the halo 
merger trees conserve mass and the mass of a branch increases 
monotonically (or stays the same) with time. As a result, the DHalo 
masses cannot be related in a simple way to other measures of mass, 
such as the number of particles identified with a friends-of-friends 
(F0F) percolation group finder. The Munich models store "cen- 
tralMvir", which is described as the "virial mass of background 
(F0F) halo containing the galaxy", and is derived from the FOF 
mass. 

The differences in these definitions of halo mass are apparent 
in Fig.Q] in which we plot the mass function of dark matter haloes 



http: //gavo .mpa-garching .mpg . de/MyMillennium/ 



http : //galaxy- catalogue . dur .ac.uk: 8888/MyMillennium/ 



3 The query used is essentially example "Gl" from the "Mainly galaxies" 
demo queries shown on the web page. 



How robust are predictions of galaxy clustering ? 3 



i mil 1 1 mil — i i i Mini — i i 1 1 him — i i inn — i i 1 1 him — r 




10 10 10 11 10 12 10 13 10 14 10 15 



M halo / h-> M 



10- 



u 



m 
C 

S— < 
CD 

a 



Bower 2006 

Font 2008 

DeLucia 2007 

Bertone 2007 




-18 -20 -22 -24 
M - 5 log(h) 



Figure 1. The mass function of dark matter haloes using the original mass 
"labels" obtained from the Millennium Archive, shown by the solid blue 
curve for the Munich models (label=centralMvir) and the solid red curve 
for the Durham models (label=mhalo ). By rescaling the Munich halo mass 
labels by A log 10 M = 0.22, as shown by the dashed blue curve, the mass 
functions agree with one another. 



using the Durham and Munich halo mass "labels". It is immedi- 
ately clear from this plot that the halo mass labels used by each 
group do not correspond to a simple particle number returned by a 
percolation group finder. The absence of a sharp cut-off in the Mu- 
nich mass function corresponding to the 20 particle limit imposed 
on the list of FOF groups stored is due to how the virial mass is 
estimated from the number of particles that the group finder says 
belong to each halo. 

At the present day the mass functions can be brought into re- 
markably close agreement with one another by rescaling the mass 
in one of the models by a constant factor. We apply this scaling to 
the Munich masses so that afterwards, haloes with the same abun- 
dance have the same masfl We need to make this rescaling as we 
plot many of our comparisons as a function of halo mass. 



3.2 Galaxy properties 

The luminosity function is the most basic description of the galaxy 
population. As such, reproducing this statistic is a primary goal 
when setting the parameters of a galaxy formation model. 

The present day luminosity functions in the r— and K-bands 
predicted by the models are plotted in Fig. [2] Note that we have 
added 5 log h to the magnitudes stored in the Millennium Archive 
for the Munich models to put them onto the same scale as the 
Durham models. The agreement between the model predictions is 
encouraging. To some extent, the model parameters have been se- 
lected to reproduce the observed luminosity function, and so one 



4 We note that a similar scaling can be performed to match the halo mass 
functions at different redshifts. However, the scaling does not work quite so 
well as it does at z = 0, and the factor required changes with redshift. 
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Figure 2. The r-band (top) and K-band (bottom) luminosity functions at z = 
0. The model predictions are shown by the different line colours and styles, 
as indicated by the key. Note that the K-band magnitude is not available in 
the Millennium Archive for the Guo et al. model, so we do not show this 
model in the lower panel. All magnitudes are on the Vega scale and include 
the effects of dust extinction. The insets show the fraction of galaxies which 
are satellites as a function of magnitude. The line styles and colours have 
the same meaning as in the main panel. 



might expect the predictions to agree even better. The models were 
not necessarily tuned to explicitly match the luminosity functions 
in these particular bands. The later versions of the Munich models 
put more emphasis on matching the inferred stellar mass function. 
Furthermore, other observables are matched at the same time as 
reproducing the luminosity function data, which may have led to 
compromises in the quality of the reproduction of the luminosity 
function. Finally, the parameter values were set by doing compar- 
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Bower et al., 2006 


9.00 


9.33 


0.070 


10.00 


9.97 


0.65 


11.00 


10.58 


7.91 


Fontetal., 2008 


9.02 


9.39 


0.090 


9.99 


9.98 


0.61 


11.04 


10.58 


7.16 


De Lucia etal., 2007 


9.17 


9.04 


0.064 


10.06 


9.57 


0.58 


11.05 


10.33 


5.90 


Guoetal., 2011 


8.90 


9.02 


0.0112 


10.06 


9.47 


0.31 


10.96 


10.24 


5.14 


Bertone et al., 2007 


9.02 


8.90 


0.006 


10.21 


9.54 


0.66 


11.07 


10.40 


8.93 



Table 1. The upper rows give the abundance of galaxies in the three "density cut" samples used in the paper. The first column below this gives the name of 
the semi-analytical model. Columns 2, 3 and 4 give the cuts applied to each model in the logarithm of stellar mass, the logarithm of the cold gas mass and the 
star formation rate, respectively for the highest density sample, density cut 1. In all cases the units of mass are /i~'M Q and the units of star formation rate are 
M©yr _1 . Columns 5-7 and 8-10 give the analogous cuts for density cuts 2 and 3 respectively. 



isons "by eye" between the model predictions and the data. It is 
not currently possible to provide a definitive list of precisely which 
datasets were used by the various authors to set the model parame- 
ters, or to specify the priorities assigned to the reproduction of dif- 
ferent datasets. This should become more transparent with future 
releases of the models, following the development of statistical ap- 
proaches to quantify goodness of fit and the weighting of datasets in 
the parameter setting proce ss dHenriques et al.ll2009t bower et alj 
l20irj ; lHenriaues et alj20"l3) . 

The inset panels in Fig. [2] show the fraction of galaxies that 
are satellites of the central galaxy within each halo as a function 
of magnitude. Again the models show similar trends, with just un- 
der half of the galaxies being satellites over most of the magnitude 
range plotted, before this value drops steadily brightwards of L,. 
Satellites and centrals are described separately in halo occupation 
distributions, so this will have implications later on. 

When we study the clustering predicted by the models, we will 
use samples defined by intrinsic properties: stellar mass, cold gas 
mass and star formation rate. The cumulative abundance of galaxies 
ranked by each of these properties in turn is shown in Fig. [3] There 
is remarkably close agreement between the distributions ranked in 
terms of stellar mass (left panel of Fig.|3}, which is surprising given 
the differences in the choice of stellar initial mass function in the 
models, which means that for a given mass of stars made, different 
amounts of light will be produced, and at least some of the models 
have been specified to reproduce observed luminosity functions. 

The model predictions agree less well when galaxies are 
ranked in terms of cold gas. The Durham models predict more 
galaxies for cold gas masses in excess of 10 9 /r'M G . This can be 
traced to less weight being given to fitting the observed gas frac- 
tions in spiral galaxies in the Bower2006 and Font2008 models. We 
note that the n ew treatment of star formation in the Durham models 
introduced bv lLagos et alj J201 1 al ibi), which distinguishes between 
atomic and molecular hydrogen in the interstellar medium, gives 
an excellent match to the observed HI mass functior0. The distri- 
butions ranked by SFR are more similar to one another than those 
for cold gas, presumably because in all cases weight was given to 
matching the observed optical colour distributions of galaxies. 

To ensure we are making a fair comparison between the mod- 
els, we define our galaxy samples by number density rather than by 
a fixed value of a property, such as the stellar mass. Hence, to ob- 

5 This model is not included in this paper because, at the time of writing, 
it was not available in the Millennium Simulation database. The model will 
be added early in 2013. 
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mass 


mass 




Density Cut 2 


Bower2006 


40% 


3% 


3% 


Font2008 


42% 


8% 


10% 


DeLucia2007 


38% 


6% 


8% 


Bertone2007 


37% 


9% 


14% 


Guo201 1 


41% 


24% 


19% 


Density Cut 3 


Bower2006 


27% 


2% 


1% 


Font2008 


26% 


6% 


2% 


DeLucia2007 


24% 


2% 


4% 


Bertone2007 


27% 


3% 


6% 


Guo201 1 


26% 


15% 


16% 



Table 2. The percentage of galaxies that are satellites in each sample. The 
upper half of the table gives the satellite percentages for density cut 2 sam- 
ples and the lower half for density cut 3. The first column gives the model 
label and columns 2, 3 and 4 give the satellite percentage for samples se- 
lected by stellar mass, cold gas mass and star formation rate, respectively. 



tain samples which contain the same number of galaxies, slightly 
different cuts on a particular property are applied to each model. 
The cuts used on each property and the number densities of the 
high, intermediate and low density samples are listed in Table 1. 
This idea of compar ing galaxy ca t alogue s at a fixed number den- 
sity was applied by iBerlind et all d2003l) and IZheng et all J2005I) 
in their HOD analysis of galaxies output by a gas dynamic simu- 
lation and a semi-analytical model. In the Berlind et al. study, the 
galaxy mass functions output by the two galaxy formation mod- 
els were very different. Nevertheless, by comparing samples as a 
fixed abundance, these authors were able to find common features 
in the HODs. The percentage of galaxies that are satellite galaxies 
is listed in Table[2] 

The finite resolution of the Millennium simulation means that 
the model predictions are incomplete below some number density. 
To investigate this we plot the results for th e Guo et al. model ob- 
tained from the Millennium II simulation teovlan-Kolchin et alj 
which has a halo mass resolution that is 125 times better 
than that of the Millennium-I run. The comparison between the 
galaxy catalogues from the Millennium-I and Millennium-II runs 
shown in Fig. [3] indicates that the Guo et al. model predictions 
from Millennium-I are robust for stellar masses and star formation 
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Figure 3. The cumulative abundance of model galaxies ranked by stellar mass (left panel), cold gas mass (middle panel) and star formation rate (right panel). 
As before, the key indicates the line colours and style used to represent each model. The three horizontal dashed lines in each panel show the three number 
densities (high, intermediate and low) used to define galaxy samples. 
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Figure 4. The ratio of stellar mass to the mass of the host dark matter halo plotted as a function of stellar mass. The halo masses from the Munich models have 
been globally rescaled as illustrated in Fig.[T] Each panel shows a different model as labelled. The horizontal blue solid line shows the maximum possible value 
of the ratio corresponding to the universal baryon fraction being converted into stellar mass in one galaxy (i.e. with no hot gas, cold gas or satellite galaxies in 
the halo). The black curves show the median (solid), 90th (dashed) and 99th (dotted) percentiles for the distribution of predicted ratios. The right-most panel 
shows the DeLucia2007 model after attempting to relabel halo masses following a postprocessing of the merger trees as explained in the text. 



rates down to our highest number density cut. The cumulative dis- 
tributions agree extremely well from the two simulations for stellar 
mass. In the case of star formation rate, the shapes of the two dis- 
tributions agree quite well, but with a small displacement. For cold 
gas, the distributions are different at low masses. For this reason, 
we omit comparisons corresponding to the first density cut in the 
case of cold gas and star formation rate, which to some extent is 
controlled by the cold gas content. Ideally this exercise should be 
repeated by comparing the Millennium-I and Millennium-II pre- 
dictions for each model in turn, as the convergence is likely to be 
model dependent. However, the Millennium-II outputs are not cur- 
rently available in the database for each model. 

Our methodology of matching the halo mass "labels" between 
models and of comparing samples with, by construction, exactly 



the same number of galaxies allows us to focus on differences in 
the way in which the models populate haloes with galaxies. 

3.3 Further adjustments to the halo masses in the Munich 
model 

We are almost in a position to compare HODs between models. 
Before we do this, we first carry out a preliminary investigation of 
how galaxy properties correlate with halo mass, by plotting the stel- 
lar mass to the host halo mass ratio in Fig. |4] This ratio is formed 
using the mass of the host halo at the present day which is the 
relevant mass for HODs (rather than the subhalo mass associated 
wit h each galaxy, whi ch is used in subhalo abundance matching 
e.g. ISimha et alJHoi3) . The median ratios of stellar mass to halo 
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Figure 5. The top row shows the HODs predicted by the models for the high (left), intermediate (middle) and low (right) density samples of galaxies ranked 
in order of descending stellar mass. Different colours correspond to different models, as indicated by the key. The DeLucia2007* model corresponds to a 
relabelling of some of the halo masses as described in the text. The satellite HOD is shown separately by a red dotted line in the case of the Bower2006 model. 
The bottom row shows the contribution to the effective galaxy clustering bias as a function of halo mass. Again the dotted red line shows the contribution to 
the bias of satellites in the Bower2006 model. Note that the DeLucia2007* model is not plotted in the lower panels as it is indistinguishable from the original 
DeLucia2007 model. 



mass are similar between the models, as shown by the solid lines in 
Fig. |U The extremes of the distribution and the outliers, are, how- 
ever, different. 

The maximum value of the ratio of stellar mass to host halo 
mass is the universal baryon fraction, since all of the models as- 
sume that dark matter haloes initially have this baryonic mass as- 
sociated with them. This extreme case corresponds to the presence 
of a single galaxy in the halo with all of the baryons in the form 
of stars, with no hot gas or cold gas. The first two panels of Fig. [4] 
show that the Durham models match this expectation with all of 
the galaxies lying below the limit indicated by the horizontal blue 
dotted line. Indeed most galaxies are far away from this line, as 
indicated by the percentile curves, which reflect s how inefficient 
galaxy form ation is at turning baryons into stars (ICole et al.ll200Tl ; 
lBaughll2006l) . 

The Munich models, on the other hand, throw up a small num- 
ber of galaxies (fewer than 1%) in which the stellar mass appears 
to exceed the universal baryon fraction, in some cases by a substan- 
tial factor, becoming comparable to the host halo mass (in fact in a 
very small number of cases the galaxy stellar masses even exceeds 
the associated halo mass). Closer inspection of the merger histories 



reveals that these galaxies, despite being labelled at the present day 
as central galaxies, reside in haloes that have been tidally stripped 
in the past. The halo spent one or more snapshots orbiting within 
a larger halo during which a substantial amount of mass was lost, 
leading to artificially high stellar mass to halo mass ratios. This sce- 
nario does not happen by construction with the Durham algorithm 
used to build merger trees (see Merson et al. 2012 and Jiang et al., 
in preparation). The stripping of mass from halos is a physical ef- 
fect. However, the extent of the mass stripping will be somewhat 
dependent on the simulation parameters. This ambiguity over the 
labelling of halo masses will have an impact on the shape of the 
predicted HOD, particularly at low halo masses. 

To enable a fair comparison between the predictions between 
the Durham and Munich models, we attempted an approximate cor- 
rection to the Munich halo masses. By examining the galaxy merger 
trees, we identified galaxies which, at the present day, are labelled 
as a central galaxy, but whose host halo was once a satellite halo 
in a more massive structure. We then trace the more massive host 
halo to the present day and assign this halo mass to the galaxy in- 
stead. The galaxy is relabelled as a satellite and is associated with 
the more massive halo, resulting in the galaxy moving down and 
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Figure 6. The predicted two-point correlation function in real-space for the stellar mass selected samples (high density - left panels; intermediate density - 
middle panels; low density - right panels). The upper panels show the correlation function and the lower panels show the correlation function multiplied by r 2 
to emphasize the differences between the predictions of the different models. Different colours and line styles represent different models, as indicated by the 
key. The dotted lines in the upper panels show f = 1, which can be taken as a measure of the correlation length ro. The dotted lines in the lower panels show 
r? for reference, to emphasize the departures from a power law. 



across in Fig. [4] We have performed this postprocessing only in the 
case of the DeLucia2007 model, which we label as DeLucia2007* 
in the right-most panel of Fig. [4] This procedure affects 3.6% of the 
galaxies in the DeLucia2007 model, altering the tail of the distribu- 
tion of the stellar mass to halo mass ratio, moving the 99% line (but 
not the 90% line) and greatly reduces the number of galaxies whose 
stellar mass exceeds the mass of baryons attached to the associated 
host halo. 



4 RESULTS 

The main results of the paper are presented in this section, follow- 
ing the preparatory work on the galaxy samples downloaded from 
the Millennium Archive, as outlined in the previous section. In Sec- 
tion 4.1, we compare the model predictions for samples defined in 
terms of stellar mass, looking at the HOD and the two-point cor- 
relation function. We examine the contribution to the correlation 
function from satellite and central galaxies, and compare the radial 
extent of galaxies in common haloes. A similar study is made for 
samples ranked by cold gas in Section 4.2 and star formation rate 
in Section 4.3. 



4.1 Stellar mass 

The HODs predicted by the models when galaxies are ranked by 
their stellar mass are plotted in the top panel of Fig. [5] The agree- 
ment between the model predictions for the high density sample 
is spectacular (left panel of Fig. |5). The Munich models display 
a slight kink at very low mean numbers of galaxies per halo com- 
pared to the Durham models. This is the regime in which only a tiny 
fraction of haloes, roughly 1 in 1000, contain a galaxy which meets 
the selection criteria. In the case of DeLucia2007*, where we have 
applied a further correction to a small fraction of halo masses as 
described in the previous section, the agreement with the Durham 
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Figure 7. The modulus of the two point correlation function predicted by 
the Bower2006 and DeLucia2007 models for the highest density stellar 
mass selected sample. The top-left panel shows the full correlation function. 
The top-right panel shows the 1-halo term, corresponding to pairs within 
the same dark matter halo. The lower panels show the 2-halo terms, with 
the correlation function of central-central pairs in the lower-left panel and 
the full 2-halo term in the lower right panel. Dotted curves show where 
the correlation function was negative before taking the modulus. For pair 
separations at which the coloured dotted curves are unity, this implies that 
if = -1 due to there being no pairs at these separations (top-right panel). 
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Figure 8. The median radius of satellites from the central galaxy in a dark 
matter halo, for the highest density sample selected by stellar mass, for the 
Bower2006 and DeLucia2007 models, as indicated by the key. The error- 
bars show the 20-80 percentile range of the distribution. The models predict 
very similar numbers of galaxies over this range of halo mass . The galaxy 
distribution has a larger radial extent in the DeLucia model. The median ra- 
dius for type 1 satellites (those which retain an identifiable subhalo; dotted 
line) and type 2 satellites (those for which the associated subhalo can no 
longer be found; dashed line) are also plotted. 



models improves and holds for all masses plotted. Here, a galaxy 
that was originally a central galaxy in a halo which has been heav- 
ily stripped is now treated as a satellite galaxy in a more massive 
halo. 

The predicted HODs are qualitatively similar for the interme- 
diate and low density samples (middle and right panels of Fig. [5}. 
However, in detail the HODs are different, particularly around the 
halo masses at which central galaxies first start to make it into the 
samples. Note that in the lowest density sample (corresponding to 
the most massive galaxies), there is no extended plateau feature 
corresponding to the situation in which all haloes contain a central 
galaxy that is included in the sample. There is a knee in the HOD 
around a halo mass of ~ 3 x 10 12 /r'M o , corresponding to the halo 
mass at which one in ten haloes contains a central galaxy which is 
sufficiently massive to be included in the sample. Satellite galaxies 
make an appreciable contribution to the HOD long before the mean 
number of galaxies per halo reaches unity, showing that the canoni- 
cal HOD model of a step function for centrals which reaches unity, 
plus a power law for satellites, is a poor description of the model 
predictions. 

The variation in the form of the predicted HOD for different 
number densities can be explained in terms of the relative impor- 
tance of AGN feedback which suppresses gas cooling in massive 
haloes. The highest number density sample (ie. corresponding to 
the lowest stellar mass cuts in each model) is dominated by galax- 
ies in haloes which are not affected by AGN feedback, since a 
mean galaxy occupation of unity is reached for haloes of mass 
~ 10 n r'M G . 

For the lower density samples (higher stellar masses), the 



form of the HOD is sensitive to the operation of AGN/radio mode 
feedback. The suppression of cooling in massive haloes changes 
the slope and scatter of th e stellar mass - halo mass relation 
dGonzalez-Perez et alj|201 lT) . With this feedback mode, the most 
massive galaxies are not necessarily in the most massive haloes. 
This is particularly true in the case of the Durham models, which 
show more scatter between stellar mass and halo mass than is dis- 
played by the Munich models (Contreras et al. in prep.). The con- 
sumption of cold gas in a starburst when a disk becomes dynami- 
cally unstable may also play a role in determining the form of the 
HOD for cold gas selected samples. Unstable disks are more com- 
mon i n the Durham models than in the Munich models jParry et al] 
l2009h . 

Differences in the HOD curves can seem quite dramatic, par- 
ticularly at low masses where the HOD makes the transition from 
zero galaxies per halo, and at high halo masses, where the HOD is 
in general a power law. However, the HOD does not give the full 
view of galaxy clustering. It is important to bear in mind that the 
number density of dark matter haloes and the clustering bias asso- 
ciated with haloes of a given mass also contribute and these quan- 
tities chan ge significantly a cross the mass range plotted in Fig. [5] 
Following lKim et alj (2003 ), we show in the bottom panel of Fig.|5] 
the contribution to the effective bias of the galaxy sample as a func- 
tion of halo mass, which takes into account the halo mass function 
and the halo bias. Note that the y-axis in the bottom panels is on a 
linear scale. The high and intermediate density samples show two 
distinct peaks in the contribution to the effective bias, with the low 
mass peak corresponding to central galaxies and the high mass peak 
to satellite galaxies. For the low density sample (bottom right panel 
of Fig. [5}, the central galaxy peak is much less distinct and satellite 
galaxies dominate the effective bias. 

The two-point correlation function of the stellar mass se- 
lected samples is plotted in Fig. [6] At larger pair separations (1 < 
(r//i~'Mpc) < 30), the correlation functions predicted by the mod- 
els are remarkably similar, as expected from the similarity in the 
effective bias plotted in Fig. [5] The variation in the predicted bias, 
averaged over pair separations of 5 - 25/7~'Mpc is less than 10%. 
At small pair separations, for the lowest density sample, the corre- 
lation functions differ by 10-30%. In the case of the intermediate 
and high density samples, the extremes of the predictions vary by a 
factor of 2-3, with more clustering predicted in the Durham models 
than in the Munich models. 

To gain further insight into the correlation functions predicted 
by the models, we examine the contributions from galaxy pairs in 
the same halo (the 1-halo term in HOD terminology) and from dif- 
ferent halos (the 2-halo term) in Fig. [7] in which we compare the 
Bower2006 and DeLucia2007 results. The top-left panel confirms 
that the correlation functions are remarkably close to one another, 
except for pair separations smaller than r ~ l/?~'Mpc. This differ- 
ence is dominated by the 1-halo term (top-right panel), which itself 
is largely determined by satellite-satellite pairs. In the Bower2006 
model the 2-halo term makes a small contribution to the amplitude 
of the correlation function at these pair separations, whereas in the 
DeLucia2007, centrals contribute no pairs at these scales. This dif- 
ference can be understood in terms of the DHalos algorithm used 
to build halo merger histories, which breaks up some FOF haloes 
into separate objects, allowing centrals to be found closer together 
than in the Munich models. 

After rescaling the halo masses and considering samples with 
the same number density of galaxies, we find that the Bower2006 
and DeLucia2007 models contain very similar numbers of satellite 
galaxies per halo, as revealed by the close agreement between the 
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HODs plotted in Fig.[5](see Table[2]). This implies that overall, the 
timescale for galaxies to merge must be similar in the models. Both 
models use analytical calculations of the time required for a satel- 
lite to merge with the central galaxy, ba s ed on the dynamic al fric- 
tion timescale (see lLacev & Cold dl993l) ; ICole et alffeOOOh ). Very 
similar expressions are used for this timescale, with an adjustable 
parameter chosen to extend the timescales in both cases to improve 
the model predictions for the bright end of the galaxy luminosity 
function. In the Munich models, the satellite orbit is followed as 
long as the associated subhalo can be resolved, then the time re- 
quired for the satellite to merge with the central is calculated ana- 
lytically. In the Durham models, the merger timescale is calculated 
as soon as the galaxy becomes a satellite, without any consideration 
as to whether or not the associated subhalo can still be identified. 

The reason for the enhanced small scale clustering in the 
Durham model must therefore be due to a difference in the spa- 
tial distribution of these satellites within dark matter haloes. Fig. [8] 
confirms this, showing the median separation between the central 
galaxy in a halo and its satellites. The median satellite radius is 
greater in the Munich model than in the Durham model. This is 
readily understood from the way in which the merger times are cal- 
culated. In the Munich model, satellite galaxies whose subhaloes 
can be resolved do not merge with the central galaxy. In the Durham 
model, some fraction of galaxies which are associated with an iden- 
tifiable subhalo will be assumed to have merged, as a result of the 
purely analytical calculation of the merger time. The spatial distri- 
bution of satellites with resolved subhaloes is more extended than 
the overall distribution of satellites. This is clear from Fig.[8]which 
shows that Type 1 satellites in the DeLucia2007 model, i.e. those 
with resolved subhaloes, have a larger median radius than satellites 
whose subhalo can no longer be identified (Type 2 satellites). There 
are around three times as many satellites with resolved subhaloes 
in the DeLucia model compared with the Bower model. 



4.2 Cold gas mass 

The models considered in this paper track the total mass in cold 
gas, combining helium and hydrogen, the latter in both its atomic 
and molecular forms. All of this material is assumed to be avail- 
able to make stars. With each episode of star formation, some cold 
gas will be ejected from the ISM and perhaps from the dark matter 
halo altogether due to supemovae. More recent models make the 
distinction between molecular and atomic hydrogen using the pres- 
sure in the mid-plane of the galactic disk and assume that only the 
molecular hydroge n is involved in star formation (see for example 
lLagosetall2011bh . 

The comparison presented in Fig. [3] between the cumulative 
cold gas mass function predicted in the Guo et al. model in the 
Millennium-I and Millennium-II simulations suggests that the pre- 
dictions have not converged for the highest density sample, corre- 
sponding to the lowest cut in cold gas mass. Hence we compare 
the predicted HODs only for the two lower density cuts in Fig. [9] 
which correspond to higher cold gas masses. There is a range of 
HOD predictions. 

The predicted HOD for central galaxies tends to show a peak 
rather than a step. This feature is particularly clear in the low- 
est density sample. This form of the HOD is due to the inclusion 
of AGN feedback in the models, which suppresses gas cooling in 
massive haloes (see the plot of cold gas mass versus halo mass in 
iKim et alj20T 2). In t he Munich models, the suppression of cooling 
comes in gradually dCroton etai1l2006l) . whereas in the Durham 
models, the suppression is total in haloes in which the AGN feed- 




Figure 9. The upper panels show the predicted HOD for the two lowest 
density samples when galaxies are ranked by their cold gas mass. Different 
line colours and styles refer to the predictions from different models, as 
indicated by the key. The lower panels show the contribution to the effective 
bias of the sample from galaxies in different mass haloes. 
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Figure 10. The predicted two-point correlation function in real-space for 
galaxy samples ranked by their cold gas mass. Different line colours and 
styles refer to the predictions from different models, as indicated by the 
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back is operative jBower et al]|2006h . We present a parametric fit 
to the form of the HOD of samples selected by cold gas mass in 
the Appendix, in w hich we advocate using the nine parameter fit of 
lGeachetalJ<2012h to describe the HOD predicted by the models. 

Satellites make up a smaller fraction of samples selected by 
cold gas mass (see Table|2j. The satellite HOD is a power-law as it 
was in the case of samples selected according to their stellar mass, 
but with a shallower slope (a < 1). There is a substantial range in 
the amplitude of the satellite HODs, which differ by a factor of 5 
between the extremes. This is consistent with the variation in the 
percentage of galaxies in each sample that are satellites, as listed 
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Figure 11. The upper panels show the predicted HOD for the two low- 
est density samples when galaxies are ranked by their star formation rate. 
Different line colours and styles refer to the predictions from different mod- 
els, as indicated by the key. The lower panels show the contribution to the 
effective bias of the sample from galaxies in different mass haloes. 



in Table [2] This is due in part to the way in which the model pa- 
rameters were set. The Bower2006 and Font2 008 models ov erpre- 
dict the observed cold gas mass function dKim etal.ll201lh . The 
Font2008 model invokes a revised model for gas cooling in satel- 
lites, in which the hot gas attached to the satellite is only partially 
stripped, with the fraction of material removed depending upon the 
orbit of the satellite. Hence, gas can still cool onto satellites in the 
Font et al. model, which explains why the satellite HOD is higher 
in this model than in Bower et af| 

The lower panels of Fig.|9]show that central galaxies dominate 
the contribution to the bias factor for cold gas selected samples. The 
peak in these contributions covers a factor of a 5 in halo mass for 
the intermediate density sample and an order of magnitude for the 
low density sample. This lack of consensus between the model pre- 
dictions is borne out in the predicted correlation functions of cold 
gas selected samples plotted in Fig. [TO] On the whole, the correla- 
tion function for cold gas samples is shallower than that for stellar 
mass samples, with a slope around —1.3 rather than -2, due to the 
less important role of satellite galaxies in the cold gas samples. 
There is a spread of a factor of a 1.3 in the clustering amplitude 
around r ~ 10/r'Mpc and of around an order of magnitude or 
more at the smallest pair separations plotted. 
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6 We note that the latest version of the Munich model jGuo et al 1 l201ll) 
includes a similar treatment of gas cooling onto satellites to that introduced 
by Font et al. However, these authors do not present a plot showing how the 
model predictions compare with cold gas data, so we cannot comment on 
the difference between the Guo2011 and Durham models for the cold gas 
mass function. 
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Figure 12. The predicted two-point correlation function in real-space for 
samples ranked by their star formation rate. Different line colours and styles 
refer to the predictions from different models, as indicated by the key. The 
lower panels show the correlation functions after multiplying by r 13 . 



4.3 Star formation rate 

Finally, we repeat the analysis of the previous section using sam- 
ples ranked by the global star formation rate in the galaxy, which 
is relevant to observational samples selected by emission line 
strength or rest-frame UV luminosity. The predicted HODs plot- 
ted in Fig. [TT] are very similar to those for the cold gas samples 
plotted in Fig. [9] In the Bower2006 and Font2008 models there is 
a direct connection between the total cold gas mass and the star 
formation rate. In the Munich models, only the gas above a surface 
density threshold is assumed to participate in star formation. The 
correlation functions plotted in Fig.[T2]also show a shallower slope 
than those obtained for the stellar mass selected samples. The pre- 
dicted clustering shows a similar spread to that displayed for cold 
gas selected samples. 



5 CONCLUSIONS 

We have tested the robustness of the predictions of semi-analytical 
models of galaxy formation for the clustering of galaxies. This 
project was made possible through the Millennium Simulation 
database, which allows the outputs of the models to be downloaded 
and analyzed by the astronomical community. 

The models tested were produced by the "Durham" and "Mu- 
nich" groups, and represent different versions of their "best bet" 
models when originally published. These codes attempt to model 
the fate of the baryonic component of the universe. The underly- 
ing physics is complex and still poorly understood despite much 
progress over the past twenty years. A lthough the two group s de- 
veloped their approaches starting from lWhite & Frenkldl99ll) , the 
implementations of the different processes which influence galaxy 
formation are quite different, as is the emphasis on which observa- 
tions should be reproduced by the models in order to set the values 
of the model paramete rs. The calculat i ons ar e set in the Millennium 
N-body simulation of Sprinael et al. (2005). However, all aspects 
of galaxy formation modelling beyond the distribution of the dark 
matter, including the construction of merger histories for halos, are 
independent. 

Given the uncertainty in the modelling of galaxy formation, 
the availability of model outputs from the same N-body simulation 
gives us an opportunity to perform a direct comparison of how the 
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models populate halos with galaxies. We have carried out this com- 
parison by looking at the model predictions for the halo occupation 
distribution (which it should be stressed is a model output) and the 
two-point correlation function. We compare samples at fixed abun- 
dances, high, medium and low, to account for any differences in the 
distributions of galaxy properties which persist between the model 
predictions. We also use a common definition of halo mass, again 
by referring to the abundance of haloes. We look at samples de- 
fined by intrinsic properties: stellar mass, cold gas mass and star 
formation rate. 

In the case of samples selected by their stellar mass, the agree- 
ment in the predicted HODs is remarkable, particularly as galaxy 
clustering is in general not one of the observables used to specify 
the model parameters. The two-halo contributions to the correla- 
tion function are the same in different models, which means that, 
under the conditions of the controlled comparison carried out here, 
the prediction for the asymptotic galaxy bias is robust. 

The one-halo term, i.e. the contribution from pairs in the 
same dark matter halo, is different. Again, under the conditions of 
our comparison, there is little difference in the number of satel- 
lite galaxies in each halo which implies that the calculation of the 
timescale required for a galaxy merger to take place is similar in the 
two sets of models. The difference lies in which galaxies are con- 
sidered for mergers. In the Durham models, any satellite can merge 
provided that the dynamical friction timescale is sufficiently short 
for the merger to have time to take place. In the Munich models, 
only those satellites which are hosted by subhaloes which can no 
longer be identified can merge. The radial distributio n of subhaloes 
is mo re extended than the distribution of dark matter dAngulo et alj 
|2009|) . so this naturally leads to larger 1-halo pair separations in the 
Munich models. It is fair to argue that here the Munich approach 
is more reasonable, though this depends on how well the subhalo 
find er works. 

iBudzvnski et al.1 d2012l) compared the radial density profile of 
galaxies in clusters using a subset of the semi-analytical models 
considered here, and argued that, for their selection, the Durham 
models produced longer merger timescales for satellites than in the 
Munich models, to explain a steeper inner profile. Our conclusion 
is different, as we argue that the merger timescales are similar but 
that the effective radial distribution of galaxies is different. This 
difference in interpretation could be due, however, to the way in 
which galaxy samples are compiled in both studies. 

The agreement between the models is less good when looking 
at samples selected by cold gas mass and star formation rate. In 
these cases the models give qualitatively similar predictions which 
differ in detail. There is a small spread in the asymptotic bias (a 
factor of ~ 1.2) and a large difference (order of magnitude) in the 
small scale correlation function. One notable difference from the 
samples selected by stellar mass is the form of the HOD predicted 
for cold gas and star formation rate samples. In the latter two cases 
the HOD for central galaxies has an asymmetrical peak, rather than 
the worn step function seen for stellar mass samples. 

The conclusion of our study is that for some galaxy proper- 
ties, different models of galaxy formation give encouragingly simi- 
lar predictions for galaxy clustering. This means that there is some 
hope of understanding galaxy bias in terms of the implications for 
galaxy formation physics, and from the point of view of removing 
it as a systematic in dark energy probes. However, this depends on 
the sample selection, and for surveys which target emission lines, 
the current predictions agree less well. Also, the absolute value of 
the bias predicted by the models may differ if an observational se- 
lection is applied to the catalogues, as would be done in a mock 



catalogue, rather than performing the controlled test we carried out 
of comparing samples at a fixed number density. It is also important 
to bear in mind that we have tested how closely the models predict 
galaxy bias in a single cosmology. In reality, we do not know the 
true underlying cosmology in the real Universe, and this will fur- 
ther complicate the attempt to extract the underlying clustering of 
the mass. As the modellin g of the gas conten t and star formation 
rate in galaxies improves jLagosetal.ll2011bl) . hopefully the pre- 
dictions for samples selected by e.g. emission line properties will 
become more robust. 
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Table Al. The best fitting parameters on applying the 5-parameter fit of 
Zheng et al. (2005) to the HOD of stellar mass selected samples predicted 
by the Bower2006 and DeLucia2007 models, for the intermediate density 
cut. 
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APPENDIX A: PARAMETRIC FITS TO THE HODS 
PREDICTED BY THE SEMI-ANALYTICAL MODELS 

Semi-analytical galaxy formation models follow the processes 
which influence the baryonic component of the Universe, in or- 
der to predict the number of galaxies hosted by dark matter haloes 
along with the properties of the galaxies. The HOD is therefore a 
natural prediction of semi-analytical models. In this appendix, we 
compare the HODs output by the semi-analytical models with ex- 
amples of the parametric forms typically used in the HOD analyses 
of galaxy clustering. Note that we fit the parametric form directly 



to the HOD output by the model and not to the HOD predictions 
for the correlation function and galaxy abu ndance, as is done when 
comparing HOD models to observations dCoorav & Shethll2002t 
IZehavi et al]2005r) . 

As we remarked in the paper, the form of the HOD predicted 
by the models is sensitive to the manner in which galaxies are se- 
lected. Following the standard practice with HOD fitting, we can 
separate the model output into contributions from central galaxies 
and satellites. The largest variations in the model predictions are 
seen in the HOD of central galaxies when different selections are 
applied. The central HOD for samples defined by stellar mass is 
almost a step function, but with a more gradual transition from a 
mean occupation of zero galaxies per halo to one galaxy per halo. 
For samples corresponding to high stellar mass galaxies, the pre- 
dicted central galaxy HOD need not reach unity, even for the most 
massive haloes. In the case of samples constructed on the basis of 
star formation rate or cold gas mass, the central HOD is peaked 
rather than a step. Such a feature would be hard to anticipate but can 
be readily understood in terms of the physics in the semi-analytical 
models. The HOD of satellite galaxies is approximately a power 
law and varies mainly in normalization between different selec- 
tions. 

The first parametrizatio n of the HOD we c onsider is the five 
parameter form proposed bv lZheng et al. I (l2005l) . The HOD is bro- 
ken up into the contributions from the mean number of satellite 
galaxies /V sat (M) and centrals N ccn (M) (i.e. the fraction of halos of 
a given mass which contain a central galaxy passing the selection 
if Af cen (M) < 1) per halo, with the overall mean number of galaxies 
per halo given by N(M) = N cen (M) + N iM (M). The central galaxy 
HOD is a rounded step function, with a gradual transition from zero 
to one central galaxy per halo on average: 



A'cen(M) = - 



1 + erf 



log M - log M mb 



(Al) 



C {ogM 

Here M is the mass of the host dark matter halo, erf(x) the error 
function, 



erf(x) : 



4= fV 



dt. 



(A2) 



M m j„ is the mass at which the half of the halos have at least one 
galaxy and o"i ogM is the width in halo mass of the transition from 
zero to one galaxy per halo. The HOD for satellite galaxies is given 
by 



iVsat(M) : 



M - M cul 
Mi 



(A3) 



where M cut is the mass below which a halo can not host a satellite 
galaxy and Mi is the mass in which a halo contains on average 1 
satellite galaxy and a is the power-law slope, which usually has a 
value close to unity. 

A schematic of the 5-parameter HOD with an explanation of 
which features the parameters control is given in the left-hand panel 
of Fig. Al. The best fitting values of the parameters of this HOD 
model are listed in Table Al for the intermediate abundance stel- 
lar mass selected samples from the Bower2006 and DeLucia2007 
models. The predicted HODs and the fits are plotted in Fig. A2. 

In the case of samples selected by cold gas mass or SFR, the 
form of the predicted HOD is very different from that for stellar 
mass sele cted samples. The HOD of central galaxies is peaked, as 
noted by iKim et alj 1201 ll) . The five parameter HOD cannot de- 
scribe this behavio ur. Instead, the nine parameter fit advocated by 
iGeach et alj j2012l) is more appropriate, with one modification. In 
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Model 


fa 




M c 


(acriogM) 


{b(T\a g M) 






<SlogM 


a 








(h- [ M e ) 








(h- [ M ) 






Cold Gas Mass 


Bower20()6 





0.675 


11.36 


0.130 


0.395 


8.00 x 10~ 3 


12.0 


0.586 


0.864 


DeLucia2007 


1.87 X 10~ 2 


0.97 


11.95 


0.34 


0.49 


3.7 x 10~ 3 


11.3 


1.126 


0.854 


SFR 


Bower2006 


6 X 1(T 4 


0.66 


11.53 


0.200 


0.380 


1.25 x 10~ 3 


11.3 


0.987 


0.942 


DeLucia2007 


3.2 X 10~ 2 


0.92 


11.86 


0.308 


0.392 


8.75 x lO" 4 


10.5 


0.205 


0.895 



Table A2. The best fitting parameters when applying the 9-parameter fit (Eqns. A3 & A4) to the HOD predicted by the Bower2()06 and DeLucia2007 models 
for intermediate density samples selected by their cold gas mass (top half) and star formation rate (bottom half). 




Figure Al. A schematic illustrating the parameteric forms used to fit the HODs predicted by the semi-analytical models. The overall HOD (blue) is separated 
into the contributions from central galaxies (red dashed) and satellites (red dotted). The left panel shows the 5-parameter fit of Zheng et al. (2005). (Eqns. 
A1-A3) and the right panel shows the 9-parameter fit of Geach et al. (2012) (Eqns. A4 and A5). The labels give the parameter names as they appear in these 
equations and indicate which features of the HOD they control. 



this case, the central HOD is now given by 
log(M/M c ) 2 



AU(M) = Ff(l-F?)exp 



Ft 



1 + erf 



2{xai ogM ) 2 
log(M/M c ) 
xc logM 



(A4) 



where F c ' are normalization factors (with values varying between 
zero and one), M c is the central value of the Gaussian. The mod- 
ification we have made is to allow the width of the Gaussian to 
be different for masses on either side of M c . The dispersion is 
xa ~\ogM - ao "io g M f° r halos with mass M < M c and ftcr logM for 
those with M > M c . This allows the best fitting central HOD to be 
an asymmetric peak. The satellite HOD is given by: 



N S AM) = F s 



1 + erf 



log(M/M raln ) 



M 



(A5) 



where F s is the mean number of satellites per halo halo at the "tran- 
sition" mass M ra j„, which corresponds roughly to the mass above 
which the mean number of galaxies per halo is dominated by satel- 
lites, S] ogM represents the width of the transition from zero satellites 
per halo to the power law, and a is the slope of the power-law which 
gives the mean number of satellites for M > M m i„ 

A schematic of the 9-parameter HOD, along with an explana- 
tion of which features the parameters control is given in the right- 
hand panel of Fig. Al. The best fitting values of the parameters of 
this HOD model are listed in Table A2 for the intermediate abun- 
dance cold gas mass and star formation rate selected samples from 
the Bower2006 and DeLucia2007 models. The predicted HODs 
and the fits are plotted in Fig. A2. 




Figure A2. The solid lines show the HOD predicted by the Bower2006 and DeLucia2007 semi-analytical models and the dashed lines show the best fitting 
parametric forms to these, with the parameters listed in Tables Al and A2. Each panel shows the predicted HODs and the fits for a different ranking of galaxies, 
based on stellar mass (left panel), cold gas mass (middle) and star formation rate (right panel). In all cases, the intermediate density sample is shown. 



